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The two-dimensional Holstein model is studied by means 
of direct Lanczos diagonalization preserving the full dynamics 
and quantum nature of phonons. We present numerical exact 
results for the single-particle spectral function, the polaronic 
quasiparticle weight, and the optical conductivity. The po- 
laron band dispersion is derived both from exact diagonaliza- 
tion of small lattices and analytic calculation of the polaron 
self-energy. 



I. INTRODUCTION 

Research on polarons has recently gained renewed in- 
terest on account of the observation of polaronic effects in 
several important classes of materials with perovskite re- 
lated structure, such as the bismuthates (Bai_ x K x Bi03), 
the high-T c cuprates (e.g., La 2 - x Sr x Cu04), the non- 
metallic nickelates (La2- x Sr x Ni04), or the colossal mag- 
netoresistive manganites (e.g., Lai_ x Ca x Mn03) [Q. De- 
spite a more than a half-century of theoretical and exper- 
imental study, our understanding of the formation and 
transport properties of polarons and bipolarons is still 
incomplete. Even the very fundamental problem of a 
single tight-binding electron coupled locally to a set of 
non-interacting Einstein oscillators (Holstein model Q|) 
has not been solved exactly, i.e., no satisfactory descrip- 
tion of the full spectral properties has been obtained so 
far. That is because the standard analytical techniques, 
based, e.g., on variational approaches || or on weak- 
coupling H and strong-coupling adiabatic and non- 
adiabatic perturbation expansions, fail to tackle this 
complicated many-body problem precisely in the physi- 
cally most important transition region, where the highly 
non-linear "self-trapping" process from an essentially 
delocalized quasi-free charge carrier at weak electron- 
phonon (EP) interaction to a "quasi-localized" small po- 
laron at strong coupling takes place. Note that the same 
problem arises in strongly coupled exciton-phonon sys- 
tems as well JjJ. 

As an attempt to close the gap between the weak and 
strong EP coupling limits, in this paper we investigate 
the Holstein model on finite lattices mainly employing 
approximation-free diagonalization techniques which, at 
the moment, provide the only reliable tool for studying 
polarons close to the crossover regime. 

Explicitly we write the Holstein Hamiltonian 
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where cj^ and are the annihilation [creation] opera- 
tors of a (spinless) fermion and a boson (phonon) at Wan- 
nier site i, respectively. In (1), the free electron transfer 
(t) is restricted to nearest-neighbour (NN) pairs ((ij)), 
the phonons are treated within harmonic approximation, 
and the EP term introduces a coupling between the local 
electron density and the dispersionsless optical phonon 
mode. Denoting by g the dimensionless EP interaction 
constant, e p — g 2 u> is the Lang-Firsov strong-coupling 
polaron binding energy with uj as the bare phonon fre- 
quency. 

The single-electron Holstein model has been studied 
extensively as a paradigmatic model for polaron for- 
mation in systems with dominant short-range electron- 
lattice interactions. The principal result is that irrespec- 
tive of the adiabaticity ratio, a — to /t, small polarons 
will be formed in a D-dimcnsional system provided that 
the two conditions A = e p /2Dt > 1 and g 2 > 1 are 
fulfilled (see Fig. 1). In the weak-coupling (adiabatic) 
limit, one expects (from scaling arguments) a quasi-free 
electron behaviour for D > 1, while in a ID system 
the carrier state becomes polaronic at arbitrary small A 
("large" polaron). Previous exact diagonalization (ED) 
work has concentrated on the ID case ]§|- |l^J^ , |l3| [L7| , 
where, however, the calculations were limited to cither 
very small clusters, rather weak EP coupling (A < 1) or 
to the adiabatic limit (w G = 0). On the other hand, the 
most interesting effects will be expected if the character- 
istic electronic and phononic energy scales are not well 
separated (A ~ 1; a ~ 1). For example, very recent nu- 
merical investigations, performed in the crossover region 
for ID chains with up to 20 sites, indicate that the po- 
laronic band structure markedly deviates from a simple 
tight -binding dispersion Jl^-[2C[| . 

Our purpose here is to extend these studies to 2D sys- 
tems and explore the spectral properties of the 2D Hol- 
stein model, whereby the focus is on the intermediate 
EP coupling and phonon frequency regime. Section 2 
presents the ED results for the single-particle spectral 
function and the optical conductivity. In Sect. 3 we 
outline an analytical approach based on an "incomplete" 
Lang-Firsov transformation, to calculate second-order 
corrections to the polaronic self-energy. Section 4 con- 
tains a discussion of the polaronic band dispersion, where 
the ED data are used to test the validity of the proposed 
theory The principal results are summarized in Sec. 5. 
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II. NUMERICAL APPROACH 

Performing direct diagonalizations of the Holstcin 
model on finite (iV-site) square lattices with periodic 
boundary conditions (PBC) , we exploit the translation^ 
invariance by the use of symmetrized basis sets in the ten- 
sorial product Hilbert space of electronic and phononic 
states 

N iKRi 

{\K- a)} = e —f]f T R, d 1 ' ' " • • ' ® il s >^» • ( 2 ) 

i—l V 

Here describes the allowed lattice translations and 

K is the total momentum of the coupled EP system. 
Since the Hilbert space associated with the phonons is 
infinite even for a finite system, we apply a trunca- 
tion procedure |2^,[l3| restricting ourselves to phononic 

states \s) p h = n*Li ~7^| i (^l) \®)ph w ith a t most M 
phonons. That means, we take into account all m— 
phonon states, where m — Yli=i n t — M, nf G [0, m] 
and 1 < s < Here ) = (M + N)\fM\N\ de- 

notes the dimension of the phononic Hilbert space. Then 
a i?-symmetrized state of the Holstein model is given as 

M S(m) 

l**> = E Y,cy\K;m,s), (3) 

m=0 s=l 

where 5(m) = (N - 1 + m)l/(N - l)!m!. 

Needless to say, that we have carefully to check for the 
convergence of both the ground-state energy,_E (M), and 
the phonon distribution, |c m | 2 (M) = J2 s \c^'^ Q \ 2 , as a 
function of the maximal number of phonons retained [fL3[ . 
Implementing an improved Lanczos algorithm on paral- 
lel computers, we are able to calculate the ground-state 
properties of systems with a total dimension of about 
10 7 . To obtain information about dynamical properties, 
we combine the Lanczos diagonalization with the Cheby- 
shev recursion and maximum entropy methods, which are 
very efficient for high-energy resolution applications [p3| . 

A. Single— particle spectral function 

By examining the dynamical properties of polarons, 
of particular importance is whether a quasiparticle-like 
excitation exists in the spectrum. In order to discuss this 
issue we have evaluated the wave-vector resolved spectral 
density function 

^)=Ei^i c ii°)i 2 ^-^) ! (4) 

n 

where c- creates an electron with momentum K; |0) is 
the vacuum state (containing no phonons). For a finite 
system, the eigenvalues E g and eigenstates \Hf g) in 



the single-particle subspace can be obtained from ED. 
As stated above, we prefer to use a moment expansion 
method in order to calculate Ag(E) directly. 

In Fig. 2, Ag(E) is displayed for two different EP cou- 
plings A = 0.125 (a-c) and A = 1.25 (d-f) (corresponding 
to the weak and intermediate-to-strong coupling situa- 
tions, respectively) at the allowed A"-vectors of the 18- 
site square lattice. To visualize the intensities (spectral 
weights) connected with the various peaks (excitations) 
we also have shown the integrated density of states 

Afg(E) = f E dE'Ag(E'). (5) 

J — oo 

Let us first consider the weak-coupling case. Owing to 
the "hybridization" of electron and phonon degrees of 
freedom, we found in each if-sector excitations being 
separated from the ground state by an energy of about 
u) . As can be seen from Afg(E), however, these predom- 
inantly phononic states have only a small admixture of 
electronic character, i.e., the main spectral weight is still 
located at energies that correspond to the bare tight- 
binding levels eg — —2t(cos K x + cosK y ). Other (satel- 
lite) peaks appear due to the existence of higher excited 
multi-phonon states. Increasing the EP interaction, we 
notice significant changes in the spectral function as e p 
exceeds half the bare bandwidth. Most notably, a band 
of polaronic states with low spectral weight evolves. This 
so-called small polaron band becomes extremely narrow 
and is well-separated from the rest of the spectra in the 
strong-coupling limit (A ^ 1). At the same time spec- 
tral weight is redistributed to the high-energy part of the 
spectra. 

To elucidate the nature of the ground state in more de- 
tail, in Fig. 3 a we have plotted the phonon distribution 
in the ground state at several EP interaction strengths 
and phonon frequencies. The results unambiguously con- 
firm the importance of multi-phonon states especially in 
the strong-coupling regime (g 2 >> 1), where |c m | 2 be- 
comes the usual Poisson distribution for small polarons. 
Here the phonons will heavily dress the electron and as 
a result the bare electron no longer behaves like a well- 
defined quasiparticle. In fact, this is exactly what has 
been observed in Fig. 2 d, where the residue of the quasi- 
particle peak at K = 0, 

2(c) = K^l4l°)lko' (6) 

is weakened. Then the question arises whether one can 
construct an appropriate quasiparticle operator having 
large spectral weight at the lowest pole in the spec- 
trum. This point was intensively discussed in the con- 
text of strongly correlated electron models (e.g., for the 
t-J model) for the case of an additional hole injected in 
the antiferromagnetic ground state of the undoped sys- 
tem p4| and, more recently, similar ideas have been ad- 
dressed for strongly coupled EP systems as well p6| , p5[ . 
Of course, in the parameter regime A > 1, a > 1 the 
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small polaron defined through the Lang-Firsov transfor- 
mation H (cf. Sec. 3) is a good quasiparticle. This was 
explicitly demonstrated in recent ED work |l6| . 

To tackle the intermediate coupling and frequency 
regime, we construct a composite polaron operator with 
momentum K on the basis of the phonon distribution 
function, 

2—1 m— 

keeping in mind that the polaronic state is characterized 
by strong on-site EP correlations. In Fig. 3 b, we have 
compared the quasiparticle weight factor Z^, computed 
by the use of the dressed electron operator (0), with the 
spectral weight simply given by (J6|). As expected, the 
results for are actually very close to those for Z^ 
at weak enough EP coupling. Clearly, in this case the 
phonon distribution is sharply peaked at the zero-phonon 
state (cf. Fig. 3 a) and the particle behaves as a nearly 
free electron. By contrast, if one increases the EP cou- 
pling strength, both electrons and phonons loose their 
own identity by forming polarons. Concomitantly we ob- 
serve less "electronic character" of the polaronic quasi- 
particle and, indeed, it can be shown analytically that (in 
the bulk limit) Z^ vanishes exponentially at very large 
couplings (A > 1). On the other hand, although care 
must be taken in the interpretation of our finite-cluster 
data for Z, the large value found for Z^ suggests that 
the proposed d-operator correctly describes the phonon 
dressing of the quasiparticle even in the transition regime. 
Accordingly the position of the maxima in the distribu- 
tion function provides an estimate of the (most probable) 
number of phonons contained in the "phonon-cloud" of 
the small polaron. 

B. Optical conductivity 

In linear response theory the real part of the optical 
conductivity is given by p4j 

Rea^) = VS(lu) + a r x % 9 {u) , (8) 

where T> denotes the Drude weight and the second term 
(sometimes called incoherent or regular part of the con- 
ductivity) can be written in a spectral representation as 

°J- W l^ E E 5(u-E m +E ). (9) 

For the Holstein model the (paramagnetic) current den- 
sity operator has the form 

jM =itY.i^+*-c\ + ^)- (10) 



According to the /-sum rule the integrated conductivity 
is related to the ground-state expectation value of the 
kinetic energy operator 

W t = -*5>t Ci + c j Ci ). (11) 

m 

By introducing the w-integrated weight S reg — S reg (oo), 

S re °(u)= f <L/<7£2V), (12) 
Jo 

and integrating in to both terms on the r.h.s. of (@), one 
arrives at 

TTP 2 T) 

S tot = — {-H t ) = -+S re <> . (13) 

As pointed out by Emin |2(| , there are two simple lim- 
its in which absorption associated with photoionization 
of Holstein polarons is well understood and the optical 
conductivity can be calculated analytically. The first 
one is the weak-coupling case, where the absorption co- 
efficient falls monotonically with increasing applied fre- 
quency. In the opposite strong-coupling (small polaron) 
limit, the optical properties are dominated by (incoher- 
ent) small polaron hopping processes accompanied by 
multi-phonon absorptions and emissions (i.e., by non- 
diagonal transitions p7j). As a result the absorption is 
peaked about u) ~ 2e p , whereby, in contrast to the case of 
large (Frohlich-type) polarons, the low-frequency side is 
evaluated above the high-frequency side of the absorp- 
tion peak p(|. Recently, for the weak-, intermediate- 
and strong-coupling regimes absorption spectra are ob- 
tained with numerical calculation of the finite-frequency 
part of the optical conductivity of (ID) finite-size Hol- 
stein models with up to four-sites pq,ttq] . 

Since it is not expected that the dimensionality (clus- 
ter size) plays a crucial role in the extreme small po- 
laron limit, in the discussion of the optical response 
we again focus on the crossover region. Results are 
presented in Fig. 4. The upper panel shows the fre- 
quency dependence of the optical conductivity (omitting 
the Drude contribution, i.e., the intra-band transitions). 
The most interesting qualitative feature of the absorp- 
tion spectrum seems to be the strongly asymmetric line- 
shape. The low-frequency peak structure may be easily 
understood in connection with the single-particle spec- 
tra. For example, an inspection of Fig. 2 d-f verifies that 
the first/second group of peaks correspond to inter-band 
transitions from the K — (0, 0) ground state to states 
with finite momenta, triggered by one/two phonon ab- 
sorption processes (u ~ E Q g — E + nuj 01 where lu q = 1.5 

and (E o R - E ) ~ 0.5 for £=(1,1), ( 2 >°)> ( 2 > 2 )' t 3 ' 1 )' 
(3,3) [see also Fig. 5 a below]). In principle such pro- 
cesses are possible even at arbitrarily small EP coupling 
because there is always a finite overlap of the correspond- 
ing wave-functions. However, the peak strength rapidly 
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decreases if the EP coupling becomes weaker. This can 
be seen from the inset of Fig. 4 a, by comparing the in- 
tegrated weight S re9 (oj) in the dissipative part of a(w). 

To analyze the role of the EP coupling strength on 
the optical response in more detail, we have displayed 
S re s together with the total sum rule (kinetic energy) in 
Fig. 4 b. S tot is directly related to the effective pola- 
ronic hopping amplitude, i p . e // = S tot /Tre 2 , which can 
be taken as a measure of a polaron's mobility p3|pl[ . 
Containing both coherent and incoherent transport pro- 
cesses, t p e ff substantially differs from the (exponen- 
tial) polaron band renormalization factors obtained an- 
alytically in the adiabatic Holstein and non-adiabatic 
Lang-Firsov cases Q (cf. Appendix). Fig. 4 b clearly 
shows the crossover from a nearly free electron, char- 
acterized by a t Pl eff that is only weakly reduced from 
its non-interacting value (ip, e //(A = 0) = 1), to a 
less mobile small (Holstein/Lang-Firsov) polaron in the 
(adiabatic/non-adiabatic) strong-coupling limit. Con- 
comitantly, in the weak-coupling regimes nearly all the 
spectral weight stays in the Drude part. As previously 
found for the (ID) four-site Holstein cluster jig ], the 
sharp decrease of S in the crossover region is driven by 
the fall of the Drude weight. By contrast the optical ab- 
sorption due to inelastic scattering processes, described 
by the regular part of the optical conductivity, becomes 
enhanced in the transition region. Of course, S reg also 
decreases as A increases in the extreme strong-coupling 
regime. 

III. THEORETICAL APPROACH 

In this section we treat the intermediate-to-strong EP 
coupling regime analytically, applying the Green's func- 
tion formalism advocated by Schnakenberg |2^] to the 
polaron problem. We start with the Holstein Hamil- 
tonian (|l|) transformed by the incomplete Lang-Firsov 
transformation U — exp{7g YliiH ~ h) c ! c i}' ^ —U^HU 

IS 

n = v E c h - E c « c h + <* + 1) . ( 14 ) 

i i 7 j i 

where 

V = -VY(2 - 7) - M , (15) 
C ii = <M,(l- 7 )(&I + & i ), (16) 
Cij = t$m = t exp{- 7 <7(&! - 6, - b) +b j )}. (17) 

In (|l4|), the (variational) parameter 7 measures the de- 
gree of the dynamical polaron effect (0 < 7 < 1); n is the 
chemical potential. 

The Green's function equations of motion deduced 
from H lead to a set of coupled equations for general- 
ized polaron (c-c) and "mixed" (c-C) Green's functions; 
the equation for the generalized polaron Green's function 



may be converted into an equation for the generalized 
self-energy and solved by iteration (see Refs. ||^ , j3l| , j32f 
for details). In the second step of iteration, the polaron 
self energy is obtained as 

Z ,(2) (miri;TO 2 T 2 ) = -(Cm 1 m 2 )(5(Ti -r 2 ) (18) 
+ E G(mT 1 ;m'T 2 ( r 2)) 

rh' in" 

(^mim'Ji^m"^}] : 

where the r-dependence of the boson operators and 
the statistical averages are to be determined using the 
Hamiltonian of independent local oscillators (the last 
term on the r.h.s. of (|l4])). Introducing the abbre- 
viations W(mim';m2m";r) for the expression in the 
square brackets of ( [l8| ) [the rhi, m', and m" refer to lat- 
tice vectors, and r = t\ — t 2 ], the Fourier transforma- 
tion of ( |l8| ) gives the self-energy equation in the space 
of Brillouin zone if -vectors and Matsubara frequencies 
ilo v = i{2v + \)tt / (3, i.e., 

Zf{iu v ) = - E (C^™,)^™™) (19) 

m2 -mi 

m 2 -mi K\s;m'm" 

1 r 

X- / dTe l{ ^^ )T W(m 1 rh';m 2 m'';T) . 
P Jo 

Owing to the properties of Cy, the first approximation 
to the self energy, given by the first term on the 

r.h.s. of (JTsj) , is non-zero only for elementary trans- 
lations h = m 2 — rhi connecting the site rhi with the 
NN sites rh 2 ; consequently, the ^-dependence of this 
term is given by 2(cosK x + cosK y ). The terms ob- 
tained in the second step of iteration yield, besides for 
the case rh 2 — rhi — 0, non-zero contributions only for 
h' = rh 2 — fh\ (vector from rhi to next NN site m 2 ) and 

rh 2 — rhi = 2h leading to the if-f unctions 4 cos K x cos K y 
and 2(cos2fC;; + cos 2K y ), respectively. In fact, denoting 

the NN sites to rhi and rh 2 by rhi + d\ and rh 2 + d%, 
respectively, the non-vanishing W{rhifh'\ rh 2 fh''; t) are 
given by 

W(m,m + di;m,m + di]T) =t 2 J 7 1 (r;2g 2 ) (20) 

for rhi — W 2 = rh and d± — d 2 , 

W(m,m + d 1 :m,rh + d 2 \r) = tHF^r; g 2 ) (21) 

for rhi = TTi2 = tu and d± 7^ d 2 , 

W(rh 1 ,m 1 +d 1 :m 2 ,rh 2 + d 2 -.T) = t*tFi{T;g 2 ) (22) 

for fh 2 — rhi = h 1 ', 2h and simultaneously rhi + di = 
rh 2 + d 2 , 
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W(m, m + d\] m, fh; r) = W(m, m; m, m + d 2 ] r) 

= r 7 (l- 7 ) £p J- 2 (r) (23) 

for fhi — fh\ = 0, 

W(mi,m 2 ;m2, m 2 ; r) = W(mi,mi;m 2 , mi; t) 

= -r 7 (l-7)£^ 2 ( r ) (24) 

for TO2 — fh\ = /i, and 

W(m, 777; 777, 777; r) = - 7) 2 w J r 3 (r) . (25) 

Here, t* — iexp{— g 2 cothi?}, g = jg and the functions 
•^1,2,3 are defined as follows 

Fl(t,k) = exp{-KCothi?} (26) 

OO 

X [2 £ I, (C) cosh[s(tf - c*t)] + / (C) - l] , 

s=l 

F 2 {t) =sinh(7?-w r)/sinli79, (27) 
F 3 (r) = f7exp{w r} + (n + 1) exp{-w D r} , (28) 

where d = (3oj /2, n = [exp{/3w } — C = K / sinhi?, 

and I s (() denote the modified Bessel functions. 

Deriving the above expressions of W, the averages (. . .) 
were evaluated applying p3| 



e Xb e' Ybt ) =exp 



-XY 



1 - exp{-/3w } 



(29) 



In particular, the functions originating from the 
^-ordered products of Cu and <Pij were calcu- 
lated as e-derivatives of the generating functional 
(T T (exp{e[6t(T 1 )+ ft^OMffa))) (at e = 0). After 
inserting (|20[)-(|25|) into (p^|), the r-integration is eas- 
ily accomplished. As a next step, the summation over 
the Matsubara frequencies iuj^ has to be carried out. For 
this aim, the Green's function in ( [l9| ) will be expressed 
by means of the spectral function A^, (w'), i.e., 



J-oo 27T K 7W s - uj' 



(30) 



and A^,(lo') = 2t7S(lu' — £gv) will be assumed, where 
= 77 + Eg), is given by the first approximation of the 
quasiparticle energy. Further approximations that will 
be made in the explicit calculations are based (i) on the 
assumption t9 3> 1 (low-temperature approximation) and 
(ii) on the limitation to negligible carrier concentration. 
In this way we obtain 



n B (suj ) + n F (£g,) 

ILUy — + SUJ 

n B {su ) + 1 - 77f(£kv) 



(31) 



s=l 



K' 



+ 



IU) V 



, - SLU 



In view of (i) and (ii), we can neglect both n B (suj ) and 
n F (t;g,) in (|3l]). The contributions to the self energy de- 
termined by the functions -7-2 (i~) and Fz{t) are treated in 
a quite analogous manner. Collecting all non-zero contri- 
butions to the r.h.s of ( |l9| ) and performing the analytical 
continuation iuj q — > Cu = uj + iS (5 — ► + ), the following 
low-temperature formula results: 



E {2) (Cu) = -2t*(cosK x + cosKy) 



K 



(32) 



s=l 



s! N ^ i<' 

K' 



~2s 



+2i*i[4 cos if^ cos + cos 2if x + cos 2K y 



00 ~o« -1 

-^i Vc (i) 



s=l 



K' 



+2**7(1 -TK^E^'^S^--' 



(1 - 7)%^ 



-4**7(1 - 7^(008^ + cos K y )— ^2 ^( w - 



£7 * (w) = I 10 — 1 denotes the first-order (Hartree- 
Fock) Green's function for independent quasiparticles 
having band energies £g; z is the NN coordination num- 
ber. 

The polaron self energy given by ( |32"| ) will be used 
to determine the polaron band dispersion. Namely, the 
quasiparticle energy is obtained as a function of K by 
solving 



u = r 1 + ReSf{uj). 



(33) 



In the strong-coupling case 7 near 1 and an extreme 
polaron band narrowing are expected; the corresponding 
(crude) approximation consists in the omission of terms 
arising owing to the "residual" EP interaction given by 
Cu and the neglect of (uj — with respect to the energy 
scale of the Poisson distribution of the oscillator energy 
fluctuations suj . Then the polaron band dispersion is 
explicitly obtained from (B2h as follows 



1 /l 



^(7)/4t=-A 7 (2-7)--(- 



- (cos K x + cos K y ) 



a \ S I K=2j 2 g 2 

2 2 

e"T 9 1 1 



(34) 



2 2a \s/k= 7 V 

x [4 cos K x cos K y 4- cos 2K X + cos 2K y ] , 

where (l/s) K means the average of s~ 1 with respect to 
the Poisson distribution with parameter n. In the limit 
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7 — > 1, (24) becomes equivalent to the result of second- 
order Rayleigh-Schrodinger perturbation theory. 

In the general case, the self-energy is a functional of 
the function fig, = u) — £g, = lu — r\ — S^j , i.e., on the 
basis of (|2|), (33) we have 



%=Re4 2) [%,]-4 ) , 
and the dispersion relation becomes 

Eg = Qg - e pl (2 - 7 ) + Sf . 



(35) 



(36) 



better agreement including higher-order corrections [g0| . 
The more important point, however, seems to be the 
use of the inc omp lete Lang-Firsov transformation in the 
derivation of (|34|). Obviously, both the Lang-Firsov for- 
mula (Eg = —£ p — 2iexp{— g 2 }(cosK x + cos K y j) and 
standard second-order strong-coupling perturbation the- 
ory (SCPT) fail to describe the band dispersion correctly. 
Moreover, as might be expected, both approximation 
considerably overestimate the renormalization of coher- 
ent bandwidth for intermediate EP couplings and phonon 
frequencies. A discussion of the strong-coupling case will 
be given separately in the Appendix. 



IV. POLARON BAND STRUCTURE 

We are now in the position to discuss the single- 
particle band dispersion of the 2D Holstein model. In 
the first place, we can extract the ^-dependence of the 
so-called "coherent" energy band, Eg, from the position 
of the lowest peak in each spectral function Ag using the 
ED data of Sec. 2.1 (cf. Fig. 2 d-f). This has be done in 
our previous work ]19fl . In the weak-coupling case and 
for phonon frequencies less than the bare electronic band- 
width, we have found a nearly unaffected tight-binding 
band near the band center and a practically flat region 
at larger momenta. As a result the coherent bandwidth 
AE = supg Eg — inf g Eg is roughly given by u> - 

The results for the intermediate-to-strong EP cou- 
pling regime are shown in Fig 5. Besides the well-known 
narrowing of the bandwidth (AE/8t — 0.0814), there are 
at least two features worth mentioning in the crossover 
region (see Fig. 5 a). Firstly, the exact band structure dif- 
fers significantly from a simple tight -binding band having 
the same bandwidth. Moreover, because of further than 
NN hopping processes generated by the EP interaction, 
the enhancement of the effective mass is weakened near 
K = 0. Secondly, it should be noted that the "flattening" 
of the band structure, discussed above for weak EP cou- 
pling, survives to relatively large EP interactions (even 
though u> > 2AE). At this point we would like to stress 
that the finite-size effects are rather small, i.e. although 
the data points belong to different system sizes we found 
a remarkably smooth behaviour of Eg. 

In order to check for the applicability of the th eory 
outlined in Sec. 3, we have compared Eg given by ( |34| ) 
with the exact results, where the parameter 7 was de- 
termined to reproduce the correct ground-state energy. 
As can be seen from Fig. 5 b, we found a surprisingly 
good agreement of the theoretical curve with the numer- 
ical data. It is natural to ask why our strong-coupling 
approach is so good even in the intermediate EP coupling 
situation. First of all, the second-order corrections to the 
self-energy involve longer-ranged polaron hopping pro- 
cesses due to (multi-) phonon absorption and emission. 
As mentioned above such processes are of particular im- 
portance in the crossover region. Thus we expect an even 



V. CONCLUSIONS 

To summarize, in this paper we have performed a 
highly complete numerical analysis of the 2D Holstein 
model studying the effect of polaron formation. The em- 
phasis was on the intermediate coupling regime, i.e., the 
transition region from nearly free electrons to small po- 
larons. The calculations are carried out by exact diago- 
nalization of finite systems including the full dynamics of 
phonons. Results are unbiased and allow to test a new 
theoretical approach applied to the calculation of the po- 
laron self-energy as well as previously proposed theories. 
Our main results are the following: 

(i) In the 2D Holstein model a continuous transition 
from a nearly free electron to a small polaron takes 
place if both criteria, A > 1 and g 2 > 1, are satis- 
fied, i.e., A and g 2 are not by themselves indepen- 
dent parameters [ p^| , ^5| . 

(ii) The small polaron state is basically a multi-phonon 
state characterized by strong on-site electron- 
phonon correlations making the particle susceptible 
to self-trapping (lj|. 

(iii) Depending on the adiabatic ratio a, there are two 
types of small polaron states in the very strong 
coupling regime: the adiabatic Holstein polaron 
and the non-adiabatic Lang-Firsov polaron (cf. 
Appendix). By the use of the correct polaron the- 
ory, i.e. the adiabatic Holstein approximation and 
the canonical Lang-Firsov approach with appropri- 
ate corrections, one obtains an excellent estimate 
of the coherent bandwidth in both adiabatic and 
non-adiabatic regimes ||ri"| , |34| 

(iv) The kinetic energy loss as a function of the EP 
coupling strength indicates that the "width" of the 
crossover region is rather small (broad) in the adi- 
abatic (non-adiabatic) regime pi|p^[ . 

(v) The single-particle spectral function shows the evo- 
lution of a well-separated narrow small polaron 
band with low spectral weight. On the other hand, 
composite particle operators, properly dressed by 



G 



dynamical phonons, give large spectral weight in 
the spectral function, i.e., they are much better fun- 
damental excitations of the systems @,|5) . 

(vi) At intermediate EP coupling strengths and phonon 
frequencies the effective polaronic band disper- 
sion deviates substantially from a simple tight- 
binding cosine band due to further than nearest- 
neighbour ranged hopping processes generated by 
the EP interaction. The proposed theoretical 
strong-coupling approach, based on an incomplete 
Lang-Firsov transformation, yields a satisfactory 
description of the polaron band structure, whereas 
the standard Lang-Firsov and strong-coupling per- 
turbation theories fail to give the correct energy 
dispersion (coherent bandwidth). 

(vi) The lineshapc of the optical absorption spectra 
is highly asymmetric in the weak-to-intermediate 
coupling regime, whereby the optical response on 
the low-energy side of the absorption peak exceeds 
that on the high-energy side p6| , |l2| . 
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APPENDIX 

Here we would like to add some comments on the na- 
ture of the small polaron state in the strong-coupling 
limit of the Holstein model. There is some confusion 
in the literature, particularly about the use of the cor- 
rect theoryof polarons (with appropriate corrections, cf. 
Refs. [p~o| p~2|] ) . This may be originated by the fact, that 
in the extreme strong-coupling limit both the adiabatic 
Holstein j|] and non-adiabatic Lang-Firsov || formu- 
lae, obtained from expansions in powers of a <C 1 [p|]35| 
and 1/A <C 1 f36| , respectively, yield the same exponen- 
tial band renormalization: exp{— g 2 }. However, the pre- 
exponential factors are different for any value of A [ fl2| , 
and the generalized Holstein formula derived by Alexan- 
drov and Mott [M also contains essential corrections of 
the exponent, g — » g 2 , at intermediate-to-strong EP 
couplings. 

These findings are supported by exact diagonaliza- 
tions, which we have performed for the ID Holstein model 
at fixed g 2 but different adiabatic ratios a. The compar- 
ison of the exact results with the theoretical approaches 
is presented in Fig. 6. In the strong-coupling anti- 
adiabatic limit (g 2 — 5, a = 10), where our theory coin- 
cides with the SCPT, we obtained an excellent agreement 
with the exact dispersion for all K. Moreover, the simple 



Lang-Firsov formula works perfectly well in the determi- 
nation of the coherent bandwidth AE^p ~ AE = 0.027. 
This is in contrast to the adiabatic strong-coupling case 
(a = 0.8, A = 2, g 2 = 5). Here the non-adiabatic Lang- 
Firsov and SCPT expressions are found to be quite in- 
adequate, e.g., they underestimate the exact bandwidth 
AE = 0.0562 by more than a factor of two. Although 
the theory of Sec. 3 improves this result to some extent 
(AE(j = 0.97) = 0.0362), a much better estimate of 
the bandwidth is achieved using the generalized Holstein 
formula j34|]: AEh = 0.0577. In any case, in real solids 
the coherence of these small-polaron band states will be 
rapidly destroyed by thermal fluctuations and/or disor- 
der effects. 
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FIG. 1. Phase diagram of the 2D Holstein model 
obtained for a 256-site square lattice using the vari- 
ational Lanczos approach proposed by the authors in 
Rcf. pH . This technique, based on an inhomogeneous 
variational Lang-Firsov transformation, correctly repro- 
duces the adiabatic (a <C 1) and anti-adiabatic (a ^S> 1), 
weak- (A <C 1) and strong-coupling (A, g 2 ^> 1) limits. 
We stress that there is no true phase transition between 
the nearly free electron and small polaron states, i.e., the 
transition line only indicates the crossover region. 

FIG. 2. Single-particle spectral function A^(E) (thin 
lines) and partial integrated density of states Afg(E) 
(bold lines) for the 2D Holstein model with 18 sites at 
A = 0.125, a = 0.8 (left column, a-c) and A = 1.25, 
a = 1.5 (right column, d-f). The if-values are given in 
units of 7r/3. Here and in the following all energies are 
measured in units of t. 

FIG. 3. Weight of the m-phonon states |c m | 2 (a) and 
wave-function renormalization Z^ [Z^] of a single elec- 
tron (b) in the K = ground state of the 2D Holstein 
model with ten sites. The residue of the "quasiparticle" 
pole shown in b was calculated using both the bare (open 
symbols) and dressed electron (filled symbols) operators. 
For further explanation see text. 

FIG. 4. Optical absorption in the 2D Holstein model. 
Using PBC, results for the regular part of the conductiv- 
ity, cr™ 9 and S reg (uj), are obtained for the 18-site lattice 
with M = 10 phonons (a), whereas the various contribu- 
tions to the sum rule are calculated on a ten-site lattice 
with M = 16 phonons (b). 

FIG. 5. Band dispersion of the 2D Holstein model 
along the highly symmetric directions of the Brillouin 
zone. In a, the exact results obtained for finite lat- 
tices with N — 16 and 18 sites are compared with a 
rescaled tight-binding band (dotted line). The solid line 
is a least-squares fit to the ED data. In b, the chain- 
dashed line gives the band structure according to the the- 
ory developed in Sec. 3. Long- and short-dashed lines 
are the results of the Lang-Firsov approximation and the 
second-order SCPT, respectively. 

FIG. 6. Small polaron band dispersion for the ID Hol- 
stein model. Exact results are compared with the predic- 
tions of different analytical approaches in the adiabatic 
(a = 0.8) and anti-adiabatic (a = 10) strong-coupling 
cases (g 2 = 5, A > 1). 
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